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© A method of measuring a thickness of a multilayered sample. 

© Spectral reflection ratios with respect to a second wavelength range are obtained and Fourier transformed to 
derive frequency converted spectrum. A power spectrum is obtained from the frequency converted spectrum to 
identified a peak which expresses interference caused by a silicon film. An approximate value d2' of the film 
thickness of the silicon film is calculated based on the peak position. After filtered by low-pass filtering, the 
frequency converted spectrum is reverse Fourier transformed to obtain spectral reflectance. From the spectral 
reflectance, theoretical spectral reflection ratios which are theoretically derived on only one transparent film of 
the thickness d3 which is formed on a silicon layer are subtracted. An approximate value d1* of the thickness of 
a silicon oxide film is then calculated from the spectral reflectance which are obtained by subtraction. 
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The present invention relates to a non-contact, non-destructive method of measuring film thicknesses of 
a multilayered sample object in which one or two transparent films are disposed on an SOI substrate which 
is formed by a body, a transparent insulation film and a monocrystalline or polycrystalline silicon film. 

In recent years, fabrication of a large scale intregration (LSI) on a Silicon/Silicon Oxide (SOI) substrate 

5 has been becoming popular. Figures 15 and 16 are cross sectional views each showing an SOI substrate 
which serves as background of the present invention. In Figures 15 and 16, a silicon oxide film (transparent 
insulation film) 1 is formed on a silicon body B and a monocrystalline silicon film 2 is formed on the silicon 
oxide film 1, thereby constituting an SOI substrate 10. As compared with a conventional LSI (which is built 
on a bulk semiconductor substrate), an LSI built on the SOI substrate 10 has better device characteristics. 

70 However, as fabrication of an LSI on an SOI substrate becomes increasingly more complex, control of 
film thicknesses needs to be more accurate than ever. In some cases, thicknesses d1 and d2 of the silicon 
oxide film 1 and the silicon film 2 which are formed on the silicon body B need to be measured. In other 
cases, where one or two transparent films are to be formed on the SOI substrate 10 during fabrication of an 
LSI. it is necessary to measure thicknesses d3 and 64 of these transparent films during the fabrication (For 

75 example, a silicon oxide film 3 needs to be formed in Figure 15 and a silicon oxide film 3 and a silicon 
nitride film 4 need to be formed in Figure 16). Further, the need is mounting to measure the thickness d1 to 
d4 at the same time. Of course, measurement of the thicknesses d1 to d4 must be non-contact and non- 
destructive since the thicknesses must be measured during fabrication of an LSI. 

Despite such needs, non-contact and non-destructive measurement of the respective thicknesses of the 

20 multilayered samples of Figures 15 and 16 is difficult. In reality, there has been no choice but to measure 
the thicknesses by destructing a multilayered sample object and observing the destructed sample object 
with an electron microscope or the like. 

To improve the situation, techniques for measuring the thickness of each layer of a multilayered sample 
object have been proposed such as that disclosed by US-A-4 999 509. According to the disclosed 

25 technique, a film thickness range of each layer is inputted in advance and the thicknesses of the respective 
layers are measured using a global optimization method and a local optimization method. 

The measurement of film thicknesses according to US-A-4 999 509, however, is not convenient for an 
operator since the operator must input a film thickness range of each layer in advance to perform the global 
optimization method. Although this disadvantage can be overcome by setting the film thickness ranges wide 

30 enough, expansion of the ranges leads to an increase in the number of computation steps, which in turn 
considerably extends the computation time. In addition, since a value calculated as a result of optimization 
largely varies depending on the starting point of the optimization (ie. the value of a film thickness of each 
layer) and other optimization parameters, the reproduction accuracy of the measurement greatly drops 
depending on setting conditions. 

as In a first aspect of the present invention there is provided a method of measuring thicknesses d1 , d2 
and d3 of a transparent insulation film, a silicon film and a transparent film, respectively, of a multilayered 
sample object which is formed by disposing the transparent film on an SOI substrate which is comprised of 
a body, the transparent insulation film and the silicon film, the method comprising: a first step of irradiating 
the multilayered sample object with light of first and second wavelength ranges to thereby measure spectral 

40 reflectance, the first wavelength range being in the ultraviolet range, the second wavelength range including 
longer wavelengths than the ultraviolet range; a second step of calculating the thickness d3 of the 
transparent film from the measured spectral reflectance of the first wavelength range; a third step of 
calculating reflectance against wave numbers at equal wave number intervals from the measured spectral 
reflectance of the second wavelength range, and frequency converting the calculated reflectance to thereby 

45 develop a frequency converted spectrum; a fourth step of obtaining a power spectrum from the frequency 
converted spectrum to find a peak which expresses interference caused by the silicon film in the power 
spectrum, and calculating an approximate value d2' of the thickness of the silicon film based on the location 
of the peak and an average refractive index of silicon in a wave number space; a fifth step of eliminating 
periodic components corresponding to a certain effective optical path or longer paths in the frequency 

so converted spectrum by low-pass filtering to thereby obtain a filtered frequency converted spectrum, then 
calculating intermediate spectral reflectance against wavelengths at equal wavelength intervals from the 
filtered frequency converted spectrum then theoretically deriving theoretical spectral reflectance of when 
only the transparent film of the thickness d3 is formed on a silicon layer, the thickness d3 being a value 
which is calculated at the second step, then subtracting the theoretical spectral reflectance from the 

55 intermediate spectral reflectance to thereby obtain final spectral reflectance for the second wavelength 
range and then calculating an approximate value dV of the thickness of the transparent insulation film from 
the final spectral reflectance; a 'sixth step of changing the thicknesses d1 and d2 from the approximate 
values d1' and d2' each by a fixed amount and calculating deviations between the theoretical spectral 
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reflectance and the measured spectral reflectance with respect to the thicknesses dl, d2 and d3 to thereby 
find a combination (dl, d2) having the minimum deviation, the thickness d3 being a value which is 
calculated at the second step; and a seventh step of performing nonlinear optimization on the thicknesses 
d1 and d2 of the combination (d1, d2) which is found at the sixth step, thereby finally determining the 
5 thicknesses d1 and d2 as the thicknesses of the transparent insulation film and the silicon film. 

In a second aspect of the present invention, preferably, the method of the first aspect further comprises 
an eighth step of increasing or decreasing the thicknesses dl and d2 each by a predetermined amount to 
tentatively determine six or more different combinations (d1, d2, d3) of the thicknesses, calculating 
deviations between the theoretical spectral reflectance and the measured spectral reflectance with respect 

70 to the six or more combinations, approximating the deviations by quadratic surface approximation utilizing 
the least square method to develop a quadratic surface function, and thereby finding and determining the 
thicknesses d1 and d2 which allow the deviation to be minimum as the thicknesses of the transparent 
insulation film and the silicon film, wherein the eighth step is repeated until the deviation between the 
theoretical spectral reflectance with respect to the thicknesses d1 and d2 calculated at the eighth step and 

15 the thickness d3 calculated at the second step and the measured spectral reflectance calculated at the first 
step becomes smaller than a certain value. 

In a third aspect of the present invention, in a method of measuring thicknesses d1, d2, d3 and d4 of a 
transparent insulation film, a silicon film and first and second transparent films, respectively, of a 
multilayered sample object which is formed by disposing the first and the second transparent films in this 

20 order on an SOI substrate which is comprised of a body, the transparent insulation film and the silicon film, 
the method comprises: a first step of irradiating said multilayered sample object with light of first and 
second wavelength ranges to thereby measure spectral reflectance, said first wavelength range being a 
ultraviolet range, said second wavelength range including longer wavelengths than the ultraviolet range; a 
second step of calculating the sum of effective optical path lengths of said first and said second transparent 

25 films from the measured spectral reflectance of said first wavelength range to find the maximum values 
d3max and d4max of the thicknesses d3 and d4, then tentatively determining the thicknesses d1 and d2 
both as zero, then changing the thicknesses d3 and d4 each by a fixed amount from zero to said maximum 
values d3max and d4max to calculate deviations between theoretical spectral reflectance and the measured 
spectral reflectance with respect to the thicknesses d1, d2, d3 and d4, then finding a combination (d3, d4) 

30 having the minimum deviation, and thereafter performing nonlinear optimization on the thicknesses d3 and 
d4 of said combination (d3, d4) and finally determining the thicknesses d3 and d4 as the thicknesses of 
said first and said second transparent films; a third step of calculating reflectance against wave numbers at 
equal wave number intervals from the measured spectral reflectance of said second wavelength range and 
frequency converting the calculated reflectance to thereby develop a frequency converted spectrum; a 

35 fourth step of obtaining a power spectrum from the frequency converted spectrum to find a peak which 
expresses interference caused by said silicon film in the power spectrum, and calculating an approximate 
value d2* of the thickness of said silicon film based on the location of the peak and an average refractive 
index of silicon in a wave number space; a fifth step of eliminating periodic components corresponding to a 
certain effective optical path or longer paths in the frequency converted spectrum by low-pass filtering to 

40 thereby obtain a filtered frequency converted spectrum, then calculating intermediate spectral reflectance 
against wavelengths at equal wavelength intervals from the filtered frequency converted spectrum, then 
theoretically deriving theoretical spectral reflectance of when only said transparent film of the thickness d3 
is formed on a silicon layer, the thickness d3 being a value which is calculated at said second step, then 
subtracting the theoretical spectral reflectance from the intermediate spectral reflectance to thereby obtain 

45 final spectral reflectance for the second wavelength range and then calculating an approximate value d1* of 
the thickness of said transparent insulation film from the final spectral reflectance; a sixth step of changing 
the thicknesses d1 and d2 from said approximate values dV and d2' each by a fixed amount, and 
calculating deviations between the theoretical spectral reflectance and the measured spectral reflectance 
with respect to the thicknesses d1, d2, d3 and d4 to thereby find a combination (dl, d2) having the 

50 minimum deviation; and a seventh step of performing nonlinear optimization on the thicknesses d1 and d2 
of said combination (dl , d2) which is found at said sixth step, thereby finally determining the thicknesses d1 
and d2 as the thicknesses of said transparent insulation film and said silicon film. 

In a fourth aspect of the present invention, preferably, the method of the third aspect further comprises 
an eighth step of tentatively determining the thicknesses d3 and d4 both as zero, then increasing or 

55 decreasing the thicknesses d3 and d4 each by a predetermined amount to tentatively determine six or more 
different combinations (d1, d2, d3, d4) of the thicknesses, then calculating deviations between the 
theoretical spectral reflectance and the measured spectral reflectance with respect to the six or more 
combinations, then approximating the deviations by quadratic surface approximation utilizing the least 
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square method to develop a quadratic surface function, and thereby finding and determining the thicknes- 
ses d3 and d4 which allow the deviation to be minimum as the thicknesses of the first and the second 
transparent films, wherein the eighth step is repeatedly performed after the second step and before the 
fourth step until the deviation between the theoretical spectral reflectance with respect to the thicknesses d1 
5 and 62 which are both zero and the thicknesses d3 and 64 having such values which are calculated at the 
eighth step and the measured spectral reflectance calculated at the first step becomes smaller than a 
certain value. 

In a fifth aspect of the present invention, method of the third aspect may alternatively further comprise a 
ninth step of increasing or decreasing the thicknesses d1 and 62 each by a predetermined amount to 

io tentatively determine six or more different combinations (d1, 62, 63) of the thicknesses, then calculating 
deviations between the theoretical spectral reflectance and the measured spectral reflectance with respect 
to the six or more combinations, then approximating the deviations by quadratic surface approximation 
utilizing the least square method to develop a quadratic surface function, and thereby finding and 
determining the thicknesses dl and 62 which allow the deviation to be minimum as the thicknesses of the 

75 transparent insulation film and the silicon film, wherein the ninth step is repeated until the deviation between 
the theoretical spectral reflectance with respect to the thicknesses d1 and 62 calculated at the ninth step 
and the thicknesses d3 and 64 calculated at the second step and the measured spectral reflectance 
calculated at the first step becomes smaller than a certain value. 

As described above, in the first aspect of the present invention, light of the first and the second 

20 wavelength ranges is irradiated upon the multilayered sample object and the spectral reflectance are 
measured. The thickness d3 of the transparent film is calculated from the measured spectral reflectance of 
the first wavelength range. Since most of light of the first wavelength range is absorbed by the silicon film 
which is disposed just under the transparent film, the spectral reflectance of the first wavelength range 
express the influence of the thickness d3 alone but not the influence of the underling layer structure 

25 including the silicon film. Hence, the thickness d3 of the transparent film is accurately calculated. 

In calculating the thicknesses d1 and 62 of the transparent insulation film and the silicon film, the 
approximate valued d1' and 62' of the thicknesses d1 and 62 are calculated by frequency analysis. This 
eliminates the conventional need to input the film thickness ranges prior to measurement of the thicknesses. 
This also reduces the number of the necessary computation steps, whereby a computation time is 

30 substantially shortened. 

Low-pass filtering is performed in calculating the approximate value dV of the thickness of the 
transparent insulation film, followed by elimination of the spectrum components of the transparent film 
having the thickness d3. Hence, even if the transparent film and the transparent insulation film are of the 
same materia! and of approximately the same thickness, the spectrum components of the transparent 

35 insulation film are accurately separated from the other spectrum components. For this reason, it is possible 
to calculate the approximate value dV of the thickness of the transparent insulation film even under such 
condition (i.e., even if the transparent film and the transparent insulation film are of the same material and of 
approximately the same thickness). 

While changing the thicknesses d1 and 62 of the transparent insulation film and the silicon film 

ao respectively from the approximate valued dV and d2 f each by a fixed amount, deviations between the 
theoretical spectral reflectance and the measured spectral reflectance with respect to the varying thicknes- 
ses d1, 62 and d3 are computed. A combination (d1, 62) which allows the deviation to be minimum is 
identified. The values d1 and 62 of the identified combination are optimized by nonlinear optimization. The 
thicknesses thus optimized are finally determined as the thicknesses d1 and 62 of the transparent insulation 

45 film and the silicon film. Thus, calculation of the thicknesses d1 and 62 is quite accurate. 

In the second aspect of the present invention, a combination (d1, 62) of the thicknesses of the 
transparent insulation film and the silicon film is found by quadratic surface approximation utilizing the least 
square method. Hence, calculation of the thicknesses d1 and 62 of the transparent insulation film and the 
silicon film is even more accurate. 

so In the third aspect of the present invention, light of the first and the second wavelength ranges is 
irradiated upon the multilayered sample object and the spectral reflectance are measured. From the 
measured spectral reflection ratios of the first wavelength range, the sum of effective optical path lengths of 
the thicknesses d3 and d4 of the first and the second transparent films is calculated, and the maximum 
values d3max and d4max of the thicknesses d3 and d4 are then found, respectively. Having tentatively 

55 determined the thicknesses dl and 62 both as zero, the thicknesses d3 and d4 are each changed by a 
fixed amount from zero in order to the maximum values d3max and d4max to calculate deviations between 
the theoretical spectral reflectance and the measured spectral reflectance with respect to the thicknesses 
d1, 62, 63 and d4. A combination (d3, d4) which allow the deviation to be minimum is then found. The 
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thickness values of the identified combination are optimized by nonlinear optimization, and the optimized 
thicknesses are finally determined as the thicknesses d3 and d4 of the first and the second transparent 
films. 

In the third aspect of the present invention, the thicknesses d1 and 62 of the transparent insulation film 
5 and the silicon film are calculated in approximately the same manner as in the first aspect of the present 
invention. 

In the fourth aspect of the present invention, after calculating the thicknesses d3 and 64 of the first and 
the second transparent films at the second step, an appropriate combination (d3, d4) of the thicknesses of 
the first and the second transparent films is identified by quadratic surface approximation utilizing the least 
70 square method. Hence, calculation of the thicknesses d3 and d4 of the first and the second transparent 
films is very accurate. 

The fifth aspect of the present invention is similar to the second aspect of the present invention in that 
an appropriate combination (d1, 62} of the thicknesses of the transparent insulation film and the silicon film 
is identified by quadratic surface approximation utilizing the least square method. Hence, calculation of the 
75 thicknesses d1 and 62 of the transparent insulation film and the silicon film is even more accurate. 

As described above, according to embodiments of the present invention, the thicknesses d1 and 62 of 
the transparent insulation film and the silicon film are calculated from the approximate values of dV and 62* 
which are calculated by frequency analysis. This eliminates a need to input the film thickness ranges prior 
to measurement of the thicknesses, and reduces the number of computation steps, thereby substantially 
20 shortening the computation time. 

Accordingly, it is an object of the present invention to offer a method of measuring thicknesses of layers 
of a multilayered sample object in which one or two transparent films are disposed on an SOI substrate 
without increasing the number of computation steps and without inputting data before making the 
measurements. 

25 In order that the present invention may be more readily understood, embodiments of the present 
invention are now described, by way of example, with reference to the accompanying drawings, in which: 
Figure 1 is a view of a film thickness measurement apparatus for use with a film thickness measurement 
method embodying the present invention; 

Figure 2 is a flow chart showing a method of measuring a film thickness of a multilayered sample object 
30 according to a preferred embodiment of the present invention; 

Figure 3 is a flow chart showing steps for measuring spectral reflection ratios of a multilayered sample 
object; 

Figure 4 is a view showing the spectral reflection ratios measured at the steps of Figure 3; 
Figure 5 is a flow chart showing computation steps for calculating the thickness d3 of a silicon oxide film; 
35 Figure 6 is a flow chart showing computation steps for calculating an approximate value d2' of the 
thickness of a silicon film from the measured reflection ratios; 
Figure 7 is a view showing an example of a power spectrum; 

Figure 8 is a flow chart showing computation steps for calculating an approximate value dV of the 
thickness of a silicon oxide film by low-pass filtering; 
40 Figure 9 is a view showing low-pass filtering; 

Figure 10 is a view showing a method of calculating the approximate value dV of the thickness of the 
silicon oxide film; 

Figure 1 1 is a flow chart showing the computation steps for finding a combination of the film thickness 
(d1, d2) which has the smallest deviation; 
45 Figure 12 is a flow chart showing computation steps for calculating the film thicknesses d1 and d2 by 
quadratic surface approximation utilizing the least square method; 

Figures 13 and 14 are flow charts showing a method of measuring a film thickness according to another 
preferred embodiment of the present invention; and 

Figures 15 and 16 are cross sectional views of examples of multilayered sample objects which serve as 
so background of the present invention. 

A.Structure of Film Thickness Measurement Apparatus 

Figure 1 is a view of a film thickness measurement apparatus for use with a film thickness measure- 
55 ment method embodying the present invention. The film thickness measurement apparatus comprises an 
illumination optical system 20 and an imaging optical system 30. The illumination optical system 20 
includes a light source 21 which is formed by a halogen lamp and a deuterium lamp so as to emit light 
within the ultraviolet range (hereinafter termed the "first wavelength range") and light of a visible range to a 
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near-infrared range (hereinafter termed the "second wavelength range") including longer wavelengths than 
the ultraviolet range. Light from the light source 21 enters the imaging optical system 30 through a 
condenser lens 22, a field stop 23 and a condenser lens 24. 

The imaging optical system 30 consists of an objective lens 31 , a beam splitter 32 and a tube lens 33. 

5 Light from the illumination optical system 20 is reflected by the beam splitter 32 and irradiated upon a 
designated illumination position IL through the objective lens 31. 

An XY stage 40 is disposed in the vicinity of the illumination position IL. A multilayered sample object 
OB such as one of those shown in Figures 15 and 16 in which the transparent films 3 and 4 are disposed 
on the SOI substrate 10 is mounted on the XY stage 40. Moving in a direction X or in a direction Y in 

70 response to a control signal received from an XY stage driving circuit (not shown), the XY stage 40 registers 
a desired surface region of the multilayered sample object OB with the illumination position IL. Although not 
shown in Figure 1 , the XY stage 40 is equipped with means for detecting the position of the XY stage 40 
(ie. an X-coordinate and a Y-coordinate). The information about the position of the XY stage 40 is supplied 
to a control unit 50 which controls the apparatus as a whole. 

75 Light reflected at the region of the multilayered sample OB (film thickness measurement region) which 
is registered with the illumination position IL is converged to a designated point on an optical axis through 
the objective lens 31 , the beam splitter 32 and the tube lens 33. A plate 62 having a pin hole 61 in the 
centre is disposed in the vicinity of the light convergent position. Of the reflected light, light passing through 
the pin hole 61 is allowed into a spectroscopic unit 70. 

20 The spectroscopic unit 70 is comprised of a concave diffraction grating 71 for separating the reflected 
light into spectrum components and an optical detector 72 for detecting spectrum components of the 
diffracted light. The optical detector 72 is formed by a photo diode array or a CCD, for example, and is 
disposed in optical conjugation with the pin hole 61. Hence, light received by the spectroscopic unit 70 is 
separated into spectrum components by the concave diffraction grating 71, and spectral signals cor- 

25 responding to the energies of relevant spectrum components are provided to the control unit 50 from the 
optical detector 72. In the control unit 50, thicknesses of a plurality of thin films which are formed in the 
multilayered sample object OB (e.g., the silicon oxide film 1, the silicon film 2, the silicon oxide film 3 and 
the silicon nitride film 4) are measured in accordance with the spectra! signals by a method which will be 
described later. The control unit 50 then outputs measurements to a display CRT 51 . 

30 As shown in Fig. 1, the control unit 50 comprises a conventional central processing unit which performs 
logic computations. Signals are transmitted between the CPU 52 and the CRT 51 and a key board 53 
through an input/output port (not shown). 

B. Operation of Film. Thickness Measurement Apparatus (Procedures of Measurement) 

35 

Fig. 2 is a flow chart showing a method of measuring a film thickness of a multilayered sample object 
according to a preferred embodiment of the present invention. A method of measuring thicknesses d1, d2 
and d3 of a multilayered sample object OB in which the silicon oxide film 3 is formed on the SOI substrate 
10 (See Fig. 15) will be hereinafter described with reference to Fig. 2. 
40 (1) At Step S1. reflection ratios Rm (X) of the multilayered sample object OB are measured on the basis 
of a calibration wafer (ie. silicon substrate) with respect to wavelengths from 200 nm to 1,000 nm. More 
specifically, measurement is performed as follows as shown in Fig. 3. First, at Step S101, the halogen 
lamp and the deuterium lamp of the light source 21 are turned on. 

Next, at Step S102, an operator sets the calibration wafer on the XY stage 40 causing light to be 
45 reflected by the calibration wafer and converged to a designated point on an optical axis through the 
objective lens 31, the beam splitter 32 and the tube lens 33. Only light passed through the pin hole 61 of 
the plate 62 is allowed into the spectroscopic unit 70 and is separated into spectrum components within 
a wavelength range from 200 nm to 1,000 nm (Step S103). At Step S104, light entering the optical 
detector 72 is photoelectrical^ converted into spectral signals which will be supplied to the control unit 
50 50 and stored as spectral data C (X) in a memory (not shown) of the control unit 50. 

Following this, the operator detaches the calibration wafer from the XY stage 40 and places the 
multilayered sample object OB which is to be measured on the XY stage 40 (Step S105). In the same 
manner as above, light reflected by the multilayered sample object OB is guided into the spectroscopic 
unit 70, is separated into spectrum components within the wavelength range from 200 nm to 1,000 nm 
55 (Step S106) and stored as spectral data M (X) in the memory of the control unit 50 (Step S107). 

Then, Steps S108 to Si 10 are repeated to read the spectral data C (X) and M (X) for every 1nm 
within the wavelength range from 200 nm to 1 ,000 nm and to calculate the reflection ratios Rm (X) from 
Eq. 1 below: 
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ton (X) = MQ1 x 100 (%) 
C(X) 



... (1) 



Data about the calculated reflection ratios (hereinafter "spectral reflection ratio data") are stored in the 
io memory (Step S109). 

Fig. 4 shows an interference waveform obtained by plotting the calculated reflection ratios Rm (X) of 
the multilayered sample object OB with respect to the wavelengths from 200nm to 1 ,000 nm against a 
wavelength X. 

(2) Proceeding to Step S2 (Fig. 2), from the spectral reflection ratios Rm(X) of the first wavelength range. 

75 the thickness d3 of the silicon oxide film 3 which is formed on the SOI substrate 10 is calculated by a 
conventional peak detection technique and the known curve fit method. Most of the light belonging to the 
first wavelength range is absorbed by the silicon film 2 which is formed just under the silicon oxide film 3 
so that the spectral reflection ratios Rm(X) of the first wavelength range express the influence of the 
thickness d3 alone but not the influence of the underling layer structure including the silicon film 2. 

20 Hence, as described below, the thickness d3 of the silicon oxide film 3 can be accurately calculated by 
the peak detection technique and the curve fit method. 

Fig. 5 is a flow chart showing the computation steps for calculating the thickness d3 of the silicon 
oxide film 3. First, at Step S201, a peak and a valley of the reflection ratios Rm (X) of the first wavelength 
range (200nm to 350nm) are found. From Eq. 2 below, an approximate value d3' of the thickness d3 of 

25 the silicon oxide film 3 is calculated. 



30 



d3 J = 0.25 x 



m-1 



nl_ 
wl 



n2_ 
w2 



35 



(2) 



where 

w1 : the wavelength of a peak (or a valley) of the shorter wavelength side 
w2: the wavelength of a peak (or a valley) of the longer wavelength side 
40 m: the total number of the peaks and the valleys in the wavelength range (i.e., w1 to w2) 

n1: the refraction index of the silicon oxide film (Si02> at the wavelength w1 
n2: the refraction index of the silicon oxide film (Si02) at the wavelength w2 
If the number of the peaks (valleys) is one or less, then d3 = 0 

At the subsequent Step S202, the thickness d1 of the silicon oxide film 1 and the thickness d2 of the 
45 silicon film 2 are set to zero. This is because most of the light belonging to the first wavelength range is 
absorbed by the silicon film 2. Upon setting d1 =d2 = 0, deviations E (200, 350 : 0, 0, d3) are calculated 
by changing the thickness d3 of the silicon oxide film 3 by a fixed amount (e.g., by 2 nm) within a certain 
range (for example, d3' ± 100 nm) from the approximate value d3\ As herein termed, "deviations E" 
refer to deviations for evaluating differences between the reflection ratios Rm (X) which are calculated 
so based on actual measurements (hereinafter "measured reflection ratios Rm (X) M ) and theoretically 
calculated reflection ratios Rc (X). Where the wavelength range is from X1 to X2 and the thicknesses of 
the silicon oxide film 1, the silicon film 2 and the silicon oxide film 3 are d1, d2 and d3, respectively, a 
deviation E at XI, X2. d1 , d2 and d3 is expressed as: 



55 
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X2 

E (XI, X2 ; dl, d2, d3) = £ (X) x (Km (X) - *c (X)F ] 

•••(3) 

In Eq. 3, W(X) is a weight function. As W(X), a weight function which is in proportion to the actual amount 
70 of received light can be used for instance. The theoretical reflection ratios Rc (X) are calculated in a 

conventional way from refraction indexes nO (X), n1 (X), n2 (X) and n3 (X) and absorption coefficients kO 

(X), k1 (X), k2 (X) and k3 (X) of each layer of the multilayered sample object OB with respect to the 

respective wavelengths, the wavelength X and the values tentatively determined as the film thicknesses 

d1, 62 and d3. This embodiment requires that the memory preliminarily stores the refraction indexes nO 
is (X), n1 (X), n2 (X) and n3 (X) and the absorption coefficients kO (X), k1 (X), k2 (X), and k3 (X) of the silicon 

body B, the silicon oxide film 1, the silicon film 2 and the silicon oxide film 3 with respect to each 

wavelength so that it is possible to read this data whenever necessary. 

From the deviations calculated in the manner above, the thickness d3 from which the minimum 

deviation is calculated is identified. Further, the deviations E are again calculated by substituting values 
20 (d3 - 2nm), d3 and (d3 + 2nm) in Eq. 4 below to develop a quadratic curve. The smallest value E of the 

resulting quadratic curve is identified and the corresponding value d3 is finally determined as the 

thickness d3 of the silicon oxide film 3. 

E = A x (d3f + B x (d3) + C (4) 

25 

where A, B and C are constants. 

(3) Next, at Step S3 (Fig. 2), the measured reflection ratios Rm (X) are frequency converted, and an 
approximate value d2' of the thickness d2 of the silicon film 2 is calculated from a peak position and a 
corresponding refraction index of the converted data. 

30 Fig. 6 is a flow chart showing computation steps for calculating the approximate value d2* of the 

thickness of the silicon film 2 from the measured reflection ratios Rm (X). In the following, the 
computation steps will be described with reference to Figs. 6 and 7. 

First, the spectral reflection ratio data with respect to the second wavelength range (e.g., 400 nm to 
1,000 nm) are read from the memory and converted into wave number data taken at equal wave number 

35 intervals (a wave number is an inverse number of a wavelength) (Step S301). This embodiment assumes 
that there are 2,048 converted data. If there is no spectral reflection ratio data which can be directly used 
as wave number data, the data which is obtained by interporating adjacent spectral reflection ratio data is 
used as desired wave number data. 

Following this, an interference waveform in the wave number space thus calculated is Fourier 

40 transformed (FFT) to obtain a Fourier transform spectrum. From the Fourier transform spectrum, a power 
spectrum as that shown in Fig. 7 is derived by plotting the absolute values of the Fourier transform 
spectrum (Step S302). The horizontal axis of the waveform of the spectrum denotes an effective optical 
path length (= refraction index x film thickness) while the vertical axis denotes the intensity of the 
spectrum components. In the graph of Fig. 7, the power spectrum includes a plurality of peaks. Because 

45 of a correlation between the refraction index and the thickness of each layer, peaks P3 and P1 related to 
the silicon oxide films 3 and 1 appear in a region in which the effective optical path is short (i.e., the left- 
hand side of Fig. 7) while a peak P2 related to the silicon film 2 appears in a region in which the effective 
optical path is long. Therefore, in this embodiment, of the sufficiently large peaks appearing for a certain 
effective optical path or longer paths (e.g., not less than 5,000 nm), a peak related to the shortest 

so effective optical path is determined as the peak P2 relating to the silicon film 2 (Step S303). 

Specifying the peak P2 as above, at Step S304, the effective optical path related to the peak P2 
( = n2 x d2) is divided by an average refractive index <n2> of the silicon film 2 for the second 
wavelength region to thereby calculate the approximate value d2' of the thickness of the silicon film 2. 

(4) Proceeding to Step S4 (Fig. 2), an approximate value dV of the thickness of the silicon oxide film 1 is 
55 calculated by low-pass filtering. More particularly, calculation follows the computation steps shown in 

Figure 8. 

First, in a manner similar to the above, the spectral reflection ratio data with respect to the second 
wavelength range (e.g., 400 nm to 1,000 nm) are converted into wave number data which are taken at 
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equal wave number intervals (Step S401), and then Fourier transformed to obtain a Fourier transform 
spectrum (Step S402). Of the Fourier transform spectrum developed at Step S402, ail data (i.e., real 
parts as well as imaginary parts) corresponding to a certain effective optical path or over (e.g., not less 
than 2,000 nm) are zeroed by low-pass filtering at Step S403, thereby eliminating all spectrum 

5 components other than components which correspond to the short effective optical paths, that is.only the 
components which relate to the silicon oxide films 1 and 3 are left. As a result, a waveform such as that 
shown in Fig. 9 is obtained (the spectrum is converted into the power spectrum for convenience sake). 
Here, it is to be noted that instead of newly developing the spectral data by Fourier transformation, the 
Fourier transform spectral data developed and stored in the memory at Step S402 can be used. 

70 Following this, at Step S404, the filtered data are reverse Fourier transformed to develop wave 

number data taken at equal wave number intervals. The data thus obtained is then converted into data 
taken at equal wavelength pitches at Step S405 to thereby find spectral reflection ratio data for the 
second wavelength region (400 nm to 1 ,000 nm). The spectral reflection ratio data thus developed does 
not include interference components caused by the silicon oxide film 2. That is, the interference 

75 waveform W13 (single dashed line of Fig. 10) is equivalent to a combination of an interference waveform 
derived on the silicon oxide film 3 (having the thickness d3) which is formed on the SOI substrate 10 and 
an interference waveform derived from the silicon oxide film 1 (having the thickness d1) which forms the 
SOI substrate 10. 

Then, at Step S406, a theoretical interference waveform WT3 (double dashed line of Fig. 10) is 
20 computed which is expected from the irradiation of the light of the second wavelength region (400nm to 

1,000nm) upon a sample which includes the silicon oxide film 3 of the thickness d3 on a silicon body. 

This computation is a conventional computation and will not be described here. 

Having developed the interference waveform W13 for the silicon oxide films 1 and 3 and the 

theoretical interference waveform WT3 for the silicon oxide 3, the theoretical interference waveform WT3 
25 is subtracted from the interference waveform W13 at Step S407. As a result, an interference waveform 

W1 (solid line of Fig. 10) expressing only the interference caused by the silicon oxide film 1 (having a 

thickness d1) which forms the SOI substrate 10 is obtained. 

Here, again by the peak detection technique, the approximate value dV of the thickness of the 

silicon oxide film 1 is calculated from the interference waveform W1 (Step S407). 
30 (5) Next, at Step S5 (Fig. 2), by changing the thicknesses d1 and d2 of the silicon oxide film 1 and the 

silicon film 2 respectively by fixed amounts (film thickness pitches) Ad1 and Ad2 from the approximate 

values dV and d2\ deviations E between the theoretical and measured spectral reflection ratios are 

calculated with respect to the varying thicknesses d1, d2 and d3. From the calculated deviations E, a 

combination of the film thicknesses (d1, d2) which allows the deviation E to be minimum. is identified. 
35 Fig. 11 is a flow chart showing computation steps for finding the combination of the thicknesses (d1, 

d2) which corresponds to the smallest deviation E. In the following the computation steps will be 

described. 

First, a deviation Emin is set at Step S501 and the film thickness pitches Ad1 and Ad2 are set at 
Step S502. Although it is desirable to take measurements from various multilayered sample objects OB 

40 and determine an appropriate value for the deviation Emin based on the measurements, the deviation 
Emin may be simply set at a sufficiently large value. With respect to the film thickness pitches Ad1 and 
Ad2, it is preferable to take a computation time into consideration in determining the film thickness 
pitches Ad1 and Ad2. In this embodiment, the film thickness pitch Ad1 is determined as 24 nm while the 
film thickness pitch Ad2 is determined as 9 nm for the reason described below. 

45 The selection of too small film thickness pitches A d1 and Ad2 will prolong the computation time. 

Therefore, the film thickness pitches Ad1 and Ad2 need to be set at reasonably large values which will 
not inconveniently extend the computation time. For instance, in an interference waveform which is 
derived for irradiation of light of 550 nm to 900 nm upon a sample which consists simply of a silicon 
oxide film (Si02) and an underlying silicon body, the total number m of peaks (or valleys) appearing in 

so the waveform is calculated as an approximate value by: 



. 900 .- 550 
m=4xnxax 

55 900 x 550 

... (5) 
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20 



where 

n: an average refractive index of the silicon oxide film with respect to light of 550 nm to 900 nm 

d: the thickness of the silicon oxide film 
As can be seen from Eq. 5, the total number m varies with the thickness d. Therefore, when the film 
thickness pitch Ad1 is large, a change in the total number m can in some cases make it impossible to 
find a minimum deviation E. To prevent this, the film thickness pitch Ad1 is determined as one-tenth of a 
change in the film thickness which is needed to increase or decrease the total number m by 1 . In other 
words, the film thickness pitch Ad1 is set as: 



Adl = 



0.1 



A , AK 900 - 550 
4 x 1.45 x 



900 x 550 



75 = 24 (nm) 



In a similar manner, the film thickness pitch Ad2 is set as: 

0.1 



Ad2 = 



4 Ann 900 - 550 
25 4 x 4.00 x 

900 x 550 

= 8.8 (nm) 



- (6) 



30 ... (7) 

Having set the deviation Emin and the film thickness pitches Ad1 and Ad2, the sequence proceeds 
to Step S503 at which the thickness d1 is changed by the film thickness pitch Ad1 in the range from (1 • 

35 a) x dV to (1 + a) x d1\ At Step S504, the thickness 62 is changed by the film thickness pitch Ad2 in 
the range from (1 - /3) x d2' to (1 + 0) x d2\ In this embodiment, the constants a and 0 are both 0.25. 

After computing a deviation E between the theoretical and measured spectral reflection ratios with 
respect to the assumed film thicknesses d1, d2 and d3, it is judged at Step S505 whether the computed 
deviation E is smaller than the deviation Emin. If "YES" at Step S505, the sequence proceeds to Step 

40 S506 and replaces the deviation Emin with the deviation E and the film thicknesses dV and d2' with the 
assumed film thicknesses d1 and d2. That is, if the newly computed deviation is smaller than the former 
deviation, the computed deviation E and the assumed film thicknesses d1 and d2 are maintained. This 
process (Step S505 and Step S506) is repeated until it is judged that a computation loops is completed 
at Step S507 and Step S508. Thus, a combination of the film thicknesses (dlmin, d2min) which causes 

45 the deviation E to be minimum is found. 

Following this, the film thickness values d1 and d2 are replaced with the film thickness values dlmin 
and d2min at Step S509. v 

(6) At Step S6 (Fig. 2), from the combination of the film thicknesses (d1 , d2), more accurate combination 
of the film thicknesses of the silicon oxide film 1 and the silicon film 2 is found by nonlinear optimization 

so (e.g., Gauss-Newton method). The Gauss-Newton method is a conventional method and will not be 
therefore described here. A wavelength range for optimization may be from 550 nm to 900 nm, for 
instance. 

(7) At Step S7 (Fig. 2), quadratic surface approximation is performed utilizing the least square method. 
From a resultant quadratic surface function, the thicknesses d1 and d2 which cause the deviation E to 

55 have a minimum value are identified. 

Fig. 12 is a flow chart showing computation steps for finding the thicknesses d1 and d2 by quadratic 
surface approximation utilizing the least square method. First, at Step S701, the thickness d1 of the 
silicon oxide film 1 and the thickness d2 of the silicon film 2 are each increased or decreased by a very 
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small amount (1 nm, for Instance) and nine combinations of the film thicknesses (d1, d2, d3) are 
tentatively determined. Deviations E between the theoretical and measured spectral reflection ratios with 
respect to the tentatively determined film thicknesses are then calculated (Table 1). 

Table 1 





dl-1 


dl 


dl+1 


d2-l 


E (XI ,X2 ; 
dl-1, d2-l, d3) 


E (XI ,X2 ; 
dl, d2-l, d3) 


E (XI ,X2 ; 
dl+1, d2-l, d3) 


d2 


E (XI ,X2 ; 
dl-1, d2, d3) 


E (XI ,X2 ; 
dl, d2, d3) 


E (XI ,X2 ; 
dl+1, d2, d3) 


d2+l 


E (XI ,X2 ; 
dl-1, d2+l, d3) 


E (XI ,X2 ; 
dl, d2+l, d3) 


E (XI ,X2 ; 
dl+1, d2+l, d3) 



At Step 702, the computed deviations E are approximated by quadratic surface approximation utilizing 
the least square method, thereby establishing a quadratic surface function as below: 

f =Ax (dl) 2 + B x (d2) 2 + C x (d1) x (d2) + D x (d1) + E x (d2) + F (8) 

where the symbols A to F are constants. 

At Step 703, from the quadratic surface function, a combination (d1, d2) which corresponds to the 
smallest deviation E is found, which is to calculate the minimum value for the deviation E. If the deviation 
E is partially differentiated by the thickness d1 or d2 and a resultant solution is zero, the deviation E has 
the minimum value E. 

ddl 



... (9) 

d 62 



... (10) 

Hence, from Eqs. 9 and 1 0, 
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2xAx (d1) + C x (d1) + 0 = 0 (11) 

2 x B x (d2) + C x (d2) + E = 0 (12) 

5 The values of the film thicknesses d1 and d2 causing the deviation E to have the minimum value are 
calculated by solving these simultaneous equations Eqs. 11 and 12. Thus, the film thicknesses are 
measured at an improved accuracy in this manner. 

(8) Next, at Step S8 (Fig. 2), the thicknesses dl, d2 and d3 calculated in the manner above are 
substituted in Eq. 3 to calculate the deviation E. It is then judged whether the deviation E thus calculated 

70 is smaller than a tolerable value. If "NO" at Step S8, the sequence returns to Step S7 to redetermine a 
combination of the film thicknesses (d1, d2). 

(9) If "YES" at Step S8, the CRT 51 displays the film thickness d1 of the silicon oxide film 1, the film 
thickness d2 of the silicon film 2 and the film thickness d3 of the silicon oxide film 3 (Step 9). This 
completes the sequence. 

75 As heretofore described, according to this embodiment, the approximate film thickness d2' of the silicon 
film 2 is calculated at Step S3 and the approximate film thickness d1' of the silicon oxide film 1 is 
calculated at Step S4, and thereafter the film thicknesses d1 and d2 are calculated from the approximate 
film thicknesses dV and d2\ Hence, it is not necessary to input the ranges of the film thicknesses before 
measuring the film thicknesses. Further, the total number of the computation steps is reduced, thereby 

20 substantially reducing the computation time. 

Although the preferred embodiment above relates only to measurement of the film thicknesses d1, d2 
and d3 of a multilayered sample object OB in which merely the silicon oxide film 3 is disposed on the SOI 
substrate 10 (Fig. 15), the present invention is also applicable to measurement of film thicknesses of a 
multilayered sample object OB in which the silicon oxide film 3 and a silicon nitride film 4 are formed on the 

25 SOI substrate 10 (see Fig. 16). Such measurement is similar to the measurement according to the 
embodiment above except for the computation steps for calculating the film thickness of the film which is 
formed on the SOI substrate 10 (namely the silicon oxide film 3 in the embodiment above). In the following, 
a method of measuring film thicknesses d1, d2, d3 and d4 of a multilayered sample object OB as that 
shown in Fig. 16 will be described. Description will be given mainly on computation steps for calculating the 

30 film thicknesses of the silicon oxide film 3 and the silicon nitride film 4 which are disposed on the SOI 
substrate 10. 

First, as in the embodiment above, spectral reflection ratios Rm (X) of a silicon substrate with respect to 
a wavelength range from 200 nm to 1 ,000 nm are measured. From spectral reflection ratios Rm (X) of the 
first wavelength range (i.e., ultraviolet range), the film thickness d3 of the silicon oxide film 3 and the film 

35 thickness d4 of silicon nitride film 4 which are formed on the SOI substrate 10 are calculated. 

Figs. 13 and 14 are flow charts showing computation steps for calculating the film thicknesses d3 and 
d4 of the multilayered sample object OB in which the silicon oxide film 3 and the silicon nitride film 4 are 
formed on the SOI substrate 10. At Step S210, a peak and a valley of the waveform, which is obtained by 
plotting the spectral reflection ratios Rm (X), in the first wavelength range (200 nm to 350 nm) are identified, 

40 and the sum od of the effective optical path lengths of the film thicknesses d3 and d4 is calculated by Eq. 
13 below: 



45 



od = 0.25 x im - 1) x wl x w2 

w2 - wl 



.. (13) 



50 Next, at Step S21 1 , the maximum value d3max for the film thickness d3 of the silicon oxide film 3 is 
calculated from Eq. 14. 



55 
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d3max = (1 + y) x od 



<n3> 

... (14) 



where 

w <n3>: an average refractive index of the silicon oxide film 3 with respect to first wavelength range 
7 : safety factor (described later) 

From Eq. 15, the maximum value d4max for the film thickness d4 of the silicon nitride film 4 is calculated. 



d4max = (1 + y) x 

<n4> 

... (15) 

20 

where 

<n4>: an average refractive index of the silicon nitride film 4 with respect to first wavelength range 
The reason why the maximum values d3max and d4max can be calculated from Eqs. 14 and 15 is because 
25 there are only a limited number of possible combinations of the film thicknesses d3 and d4 once the sum 
od of the effective optical path lengths of the film thicknesses d3 and d4 is known. In other words, a 
combination (d3, d4) is determined by finding the film thicknesses d3 and d4 which satisfy: 

(1 - y ) x od Z (<n3> x d3 + <n4> x d4) S (1 + 7 ) x od (16) 

30 

where 7, the safety factor regarding differences between measurements and calculated values (due to the 
NA and the refractive index of the objective lens, noise, etc.), is equal to or larger than 0 but smaller than 1 . 
In this embodiment, 7 = 0.2. 

At Step S212, film thickness pitches Ad3 and Ad4 are determined. The film thickness pitches Ad3 and 

35 Ad4 are set at the most appropriate values in the manner as that in Step S502. 

After setting the deviation Emin (Step S213), the film thickness d3 is changed by the pitch Ad3 within a 
range from zero to the value d3max at Step S21 4 while the film thickness d4 is changed by the pitch Ad4 
within a range from zero to the value d4max at Step S215. During this, it is judged at Step S216 whether 
Inequality 16 is satisfied. At Step 217, deviations E between the theoretical and the measured spectral 

40 reflection ratios are computed with respect to the film thicknesses d1 ( = 0), d2 ( = 0), d3 and d4, and it is 
judged whether the computed deviations E are each smaller than the deviation Emin. 

If "YES" at Step S217, the sequence proceeds to Step S218 at which the deviation Emin is replaced 
with the computed deviation E and the film thicknesses d3min and d4min are replaced with the assumed 
film thicknesses d3 and d4. This process (Step S216, Step S217 and Step S218) is repeated until it is 

45 judged at Step S219 and Step S220 that a computation loop is completed. A combination of the film 
thicknesses (d3min, d4min) which causes the deviation E to be minimum is found in this manner. 

The film thicknesses d3min and d4min are then replaced with the film thicknesses d3 and d4. From the 
combination of the film thicknesses (d3, d4) thus found, more accurate combination of the film thicknesses 
(d3, d4) is found by nonlinear optimization (e.g., Gauss-Newton method) (Step S221). 

50 At Steps S222 to S224, quadratic surface approximation is performed utilizing the least square method. 
From the resultant quadratic surface function, the thicknesses d3 and d4 which allow the deviation £ 
between the theoretical and measured spectral reflection ratios to be minimum are identified. More 
specifically, the film thickness d3 of the silicon oxide film 3 and the film thickness d4 of the silicon nitride 
film 4 are each increased or decreased by a very small amount (2nm, for example) at Step S222 so that 

55 nine different combinations (0, 0, d3, d4) are determined with respect to the thicknesses d3 and d4. 
Deviations E between the theoretical and measured spectral reflection ratios are then calculated for the 
respective combinations (Table 2). 
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Table 2 



70 



75 



20 





d4-2 


d4 


d4+2 


d3-2 


E(X1 ,X2 ; 
0, 0, d3-2, d4-2) 


E (XI ,X2 ; 
0, 0, d3-2, d4) 


E (XI ,X2 ; 
0, 0, d3-2, d4+2) 


d3 


E (XI ,X2 ; 
0, 0, d3, d4-2) 


E (XI ,X2 ; 
0, 0, d3, d4) 


F (X\ X2 ; 
0, 0, d3, d4+2) 


d3+2 


E (XI ,X2 ; 
0, 0, d3+2, d4-2) 


E (XI ,X2 ; 
0, 0, d3+2, d4) 


E (XI ,X2 ; 
0, 0, d3+2, d4+2) 
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thereof. 
Claims 

5 1. A method of measuring thicknesses d1. d2 and d3 of a transparent insulation film (1), a silicon film (2) 
and a transparent film (3), respectively, of a multilayered sample object which is formed by disposing 
said transparent film (3) on a substrate (10) which is comprised of a body (B), said transparent 
insulation film (1) and said silicon film (2), said method comprising: 

a first step of irradiating said multilayered sample object with light of first and second wavelength 

70 ranges to measure spectral reflectance, said first wavelength range being within the ultraviolet range, 
said second wavelength range including longer wavelengths than the ultraviolet range; 

a second step of calculating the thickness d3 of said transparent film (3) from the measured 
spectral reflectance of the first wavelength range; 

a third step of calculating reflectance against wave numbers at equal wave number intervals from 

75 the measured spectral reflectance of the second wavelength range, and frequency converting the 
calculated reflectance to develop a frequency converted spectrum; 

a fourth step of obtaining a power spectrum from the frequency converted spectrum to find a peak 
which expresses interference caused by said silicon film (2) in the power spectrum, and calculating an 
approximate value d2* of the thickness of said silicon film (2) based on the location of the peak and an 

20 average refractive index of silicon in a wave number space; 

a fifth step of eliminating periodic components corresponding to a certain effective optical path or 
longer paths in the frequency converted spectrum by low-pass filtering to obtain a filtered frequency 
converted spectrum, then calculating intermediate spectral reflectance against wavelengths at equal 
wavelength intervals from the filtered frequency converted spectrum, then theoretically deriving theo- 

25 retical spectral reflectance when only said transparent film (3) of the thickness d3 is formed on a silicon 
layer, the thickness d3 being a value which is calculated at said second step, then subtracting the 
theoretical spectral reflectance from the intermediate spectral reflectance to obtain a final spectral 
reflectance for the second wavelength range and then calculating an approximate value dV of the 
thickness of said transparent insulation film (1 ) from the final spectral reflectance; 

30 a sixth step of changing the thicknesses d1 and d2 from said approximate values dV and d2' each 

by a fixed amount and calculating deviations between the theoretical spectral reflectance and the 
measured spectral reflectance with respect to the thicknesses d1, d2 and d3 to find a combination 
(d1,d2) having the minimum deviation, the thickness d3 being a value which is calculated at said 
second step; and 

35 a seventh step of performing nonlinear optimization on the thicknesses d1 and d2 of said 

combination (d1. d2) which is found at said sixth step, thereby determining the thicknesses d1 and d2 
as the thicknesses of said transparent insulation film (1 ) and said silicon film (2). 

2. A method according to claim 1, further comprising an eighth step of increasing or decreasing the 
40 thicknesses d1 and d2 each by a predetermined amount to determine tentatively six or more different 

combinations (d1, d2, d3) of the thicknesses, calculating deviations between the theoretical spectral 
reflectance and the measured spectral reflectance with respect to said six or more different combina- 
tions, approximating the deviations by quadratic surface approximation utilizing the least square method 
to develop a quadratic surface function, and thereby finding and determining the thicknesses dl and d2 
45 which allow the deviation to be minimum as the thicknesses of said transparent insulation film (1) and 
said silicon film (2), wherein said eighth step is repeated until the deviation between the theoretical 
spectral reflectance with respect to the thicknesses dl and d2 calculated at said eighth step and the 
thickness d3 calculated at said second step and the measured spectral reflectance calculated at said 
first step becomes smaller than a certain value. 

50 

3. A method of measuring thicknesses dt, d2, d3 and d4 of a transparent insulation film (1), a silicon film 
(2) and first and second transparent films (3,4), respectively, of a multilayered sample object which is 
formed by disposing said first and said second transparent films (3,4) in this order on substrate (10) 
which is comprised of a body (3), said transparent insulation film (1) and said silicon film (2), said 

55 method comprising: 

a first step of irradiating said multilayered sample object with light of first and second wavelength 
ranges to measure spectral reflectance, said first wavelength range being within the ultraviolet range, 
said second wavelength range including longer wavelengths than the ultraviolet range; 
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a second step of calculating the sum of effective optical path lengths of said first and second 
transparent films (3,4) from the measured spectral reflectance of said first wavelength range to find the 
maximum values d3max and d4max of the thickness d3 and d4, then tentatively determining the 
thicknesses d1 and 62 both as zero, then changing the thicknesses d3 and d4 each by a fixed amount 
from zero to said maximum values d3max and d4max to calculate deviations between theoretical 
spectral reflectance and the measured spectral reflectance with respect to the thicknesses d1, 62, 63 
and d4, then finding a combination (d3,d4) having the minium deviation, and thereafter performing 
nonlinear optimization on the thicknesses d3 and d4 of said combination (d3, d4) and finally determin- 
ing the thicknesses d3 and d4 as the thicknesses of said first and said second transparent films (3,4); 

a third step of calculating reflectance against wave numbers at equal wave number intervals from 
the measured spectral reflectance of said second wavelength range and frequency converting the 
calculated reflectance to thereby develop a frequency converted spectrum; 

a fourth step of obtaining a power spectrum from the frequency converted spectrum to find a peak 
which expresses interference caused by said silicon film (2) in the power spectrum, and calculating an 
approximate value d2' of the thickness of said silicon film (2) based on the location of the peak and an 
average refractive index of silicon in a wave number space; 

a fifth step of eliminating periodic components corresponding to a certain effective optical path or 
longer paths in the frequency converted spectrum by low-pass filtering to obtain a filtered frequency 
converted spectrum, then calculating intermediate spectral reflectance against wavelengths at equal 
wavelength intervals from the filtered frequency converted spectrum, then theoretically deriving theo- 
retical spectral reflectance of when only said transparent film (3) of the thickness d3 is formed on a 
silicon layer (2), the thickness d3 being a value which is calculated at said second step, then 
subtracting the theoretical spectral reflectance from the intermediate spectral reflectance to thereby 
obtain final spectral reflectance for the second wavelength range and then calculating an approximate 
value dV of the thickness of said transparent insulation film (1) from the final spectral reflectance; 

a sixth step of changing the thicknesses dl and d2 from said approximate values dV and d2' each 
by a fixed amount, and calculating deviations between the theoretical spectral reflectance and the 
measured spectral reflectance with respect to the thicknesses d1, d2, d3 and d4 to thereby find a 
combination (d1, d2) having the minimum deviation; and 

a seventh step of performing nonlinear optimization on the thicknesses d1 and d2 of said 
combination (d1, d2) which is found at said sixth step, thereby finally determining the thicknesses d1 
and d2 as the thicknesses of said transparent insulation film (1) and said silicon film (2) 

A method according to Claim 3, further comprising an eighth step of tentatively determining the 
thicknesses d3 and d4 both as zero, then increasing or decreasing the thicknesses d3 and d4 each by 
a predetermined amount to tentatively determine six or more different combinations (d1, d2, d3, d4) of 
the thicknesses, then calculating deviations between the theoretical spectral reflectance and the 
measured spectral reflectance with respect to said six or more combinations, then approximating the 
deviations by quadratic surface approximation utilizing the least square method to develop a quadratic 
surface function, and thereby finding and determining the thicknesses d3 and d4 which allow the 
deviation to be minimum as the thicknesses of said first and said second* transparent films (3,4), 
wherein said eighth step is repeatedly performed after said second step and before said fourth step 
until the deviation between the theoretical spectral reflectance with respect to the thicknesses d1 and 
d2 which are both zero and the thicknesses d3 and d4 having such values which are calculated at said 
eighth step and the measured spectral reflectance calculated at said first step becomes smaller than a 
certain value. 

A method according to claims 3 or 4 comprising a further step of increasing or decreasing the 
thicknesses d1 and d2 each by a predetermined amount to tentatively determine six or more different 
combinations (d1, d2, d3) of the thicknesses, then calculating deviations between the theoretical 
spectral reflectance and the measured spectral reflectance with respect to said six or more combina- 
tions, then approximating the deviations by quadratic surface approximation utilizing the least square 
method to develop a quadratic surface function, and thereby finding and determining the thicknesses 
d1 and d2 which allow the deviation to be minimum as the thicknesses of said transparent insulation 
wherein said further step is repeated until the deviation between the theoretical spectral reflectance with 
respect to the thicknesses d1 and d2 calculated at said ninth step and the thicknesses d3 and d4 
calculated at said second step and the measured spectral reflectance calculated at said first step 
becomes smaller than a certain value. 
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A method according to any preceding claim, wherein said nonlinear optimization is a Gauss-Newton 
method. 

A method according to any preceding claim, wherein the substrate (10) is an SOI substrate 



17 



EP 0 639 753 A2 



FIG. 1 




18 



EP 0 639 753 A2 



F I G. 2 



( START ^ 



MEASURE SPECTRAL REFLECTION RATIOS 



CALCULATE d3 BY PEAK DETECTION 
AND CURVE FITTING 



-S 1 



~S 2 



CALCULATE APPROXIMATE d2 BY 
FREQUENCY CONVERSION 



CALCULATE APPROXIMATE dl BY 
LOW-PASS FILTERING 



-S3 



~S4 



CALCULATE dl & d2 HAVING MINIMUM 
DEVIATION E IN RANGES (l±a)dl & 
(\±/3 ) d2 



OPTIMIZE dl & d2 BY .OPTIMIZATION 



APPROXIMATE dl & d2 BY QUADRATIC 
SURFACE APPROXIMATION 




E < TOLERANCE 
YES 



DISPLAY FINAL dl,d2,d3 



-S 5 



-S6 



-S 7 



■S9 



± 

( END ) 



Patent provided by Sughrue Mion, PLLC - http://www.sughrue.com 
19 



EP 0 639 753 A2 



FIG. 3 



Q START ^ 



HALOGEN LAMP AND 
DEUTERIUM LAMP ON 



SET CALIBRATION WAFER 



SI 01 
S102 



SEPARATE REFLECTED LIGHT INTO 
SPECTRUM COMPONENTS (200 -KOOOnm) 



-S103 



STORE PH0T0ELECTR1CALLY CONVERTED 
SPECTRAL DATA C(A) IN MEMORY 



- SI 04 



SET MULT I LAYERED SAMPLE OBJECT 



SEPARATE REFLECTED LIGHT INTO 
•SPECTRUM COMPONENTS (200 - KOOOnm) 



- SI 05 
-S106 



STORE PHOTOELECTR I CALLY CONVERTED 
SPECTRAL DATA M ( A ) IN MEMORY |~ S 1 0 7 



CHANGE a BY ium PITCH FROM 
200 - KOOOnm 



S108 



v 


CALCULATE 


AND 


STORE 


SPECTRAL 


REFLECT 


ION 


RATIO 


DATA 




- SI 09 



20 



EP 0 639 753 A2 



FIG. 4 




21 



EP 0 639 753 A2 



FIG. 5 



Q START ) 

1 

FIND PEAK AND VALLEY FOR Rm ( A ) 
AT 200 - 350nm, 

CALCULATE APPROXIMATE THICKNESS 
d3' OF SILICON OXIDE FILE 3 



1 

CHANGE d3 BY 2nm IN RANGE 
d3' ± lOOnm AND FIND d3 HAVING 

MINIMUM DEVIATION E (200. 350; 
0, 0, d3) 



DEVELOP QUADRATIC CURVE FOR 
THREE DEVIATIONS 
E (200,350 ; 0,O,d3-2) 
E (200,350 ; 0,0, d3) 
E (200.350 ; 0.0,d3+2) 
DETERMINE SMALLEST d3 IN 
THE CURVE AS FINAL d3 



-S2 0 3 



± 

( RETURN ^ 



22 



EP 0 639 753 A2 



FIG. 6 



^ START ^ 



CONVERT SPECTRAL REFLECTION RATIO 
DATA FOR 400 - l.OOOnm INTO 2,408 
EQUAL WAVE NUMBER DATA 



-S3 0 1 



TRANSFORM 2,408 EQUAL WAVE NUMBER 
DATA BY FFT INTO POWER SPECTRUN 



-S302 



SPECIFY PEAK OF CERTAIN VALUE OR 
MORE APPEARING FOR 5,000ni OR LONGER 
EFFECTIVE OPTICAL PATH LENGTHS 
AS PEAK P2 



-S3 0 3 



CALCULATE APPROXIMATE VALUE 42, 
OF THICKNESS 41 FROM PEAK POSITION 
AND AVERAGE REFRACTUVE INDEX OF Si 
IN WAVE NUMBER REGION 



-S304 



Q RETURN J 



23 



EP 0 639 753 A2 



FIG. 7 




- o.o 

0.0 



50.0 

15.0(i?5000nm) (x 333.3 

EFFECTIVE OPTICAL PATH LENGTH 



Patent provided by Sughrue Mion, PLLC - http://www.sughrue.com 
24 



EP 0 639 753 A2 



FIG. 8 



START ^ 



CONVERT SPECTRAL REFLECTION RATIO 
DATA FOR 400 - l.OOOnm INTO 2.408 
EQUAL WAVE NUMBER DATA 



-S4 0 1 



TRANSFORM 2,408 EQUAL WAVE NUMBER 
DATA BY FFT INTO POWER SPECTRUM 



~S4 0 2 



DETERMINE ALL DATA FOR 2,000nm OR 
LONGER EFFECTIVE OPTICAL PATH 
LENGTHS AS ZERO 



~S4 0 3 



REVERSE FFT 



-S4 0 4 



CONVERT 2,408 EQUAL WAVE NUMBER DATA 
INTO SPECTRAL REFLECTION RATIO DATA 
OF 400 - KOOOnm -+ WAVEFORM (W13) 



-S4 0 5 



DEVELOP THEORETICAL INTERFERENCE 
WAVEFORM FOR 400 - l,000nm — 
WAVEFORM (WT3) 



~S4 0 6 



SUBTRACT WAVEFORM (W 13) FROM WAVEFORM 
(WT3) TO DERIVE WAVEFORM Wl, DETECT 
PEAK IN WAVEFORM Wl, CALCULATE 
APPROXIMATE VALUE dT OF THICKNESS dl 



-S4 0 7 



Q RETURN ^ 



25 



t 



EP 0 639 753 A2 



FIG. 9 





20.0 r 


if) 




\— 


18.0 - 


LU 




O 


16.0- 


CL 




O 


14.0 - 


o 






12.0- 


ZD 


q: 


10.0- 


i — 


o 


LU 




CL 
CO 


8.0- 


LL 

o 


6.0- 


>- 




h- 


4.0- 


CO 




LU 


2.0- 


I— 




2: 


0.0- 




20.0 30.0 40.0 50.0 

(X 333.3 nm) 

EFFECTIVE OPTICAL PATH LENGTH 



26 



EP 0 639 753 A2 



FIG. 10 




27 



EP 0 639 753 A2 



FIG. 11 



f START ^ 



DETERMINE Emin 



-S50 1 



DETERMINE FILM 
THICKNESS PITCHES 

Adl^24nm 
A d 2 — 9 n m 



-S 5 0 2 



CHANGE dl BY Adl FROM 
(l-a)dl' TO (l+a)dl' 



-S 503 



CHANGE d2 BY Ad2 FROM 
(l-/?)d2' TO (1+ >S )d2* t~ S504 



S5 0 5 



Emin«-£ 
d lmiiv-d 1 
d 2min«-d 2 




V 



S506 



S508 







d 1—6 lmin 
d 2«-d 2min 


1 


f 



~S 5 0 9 



Q RETURN ^ 



28 



EP 0 639 753 A2 



FIG. 12 



Q START ^ 



CHANGE dl & d: 
DETERMINE DEVI A 
COMB IN/ 


I BY ±lnn TO 
TIONS FOR NINE 
\TIONS 


> 


t 


QUADRATIC SURFACE APPROXIMATION 
USING LEAST SQUARE METHOD ON 
NINE DATA 




r 


FIND dl & d2 HAVING SMALLEST 
QUADRATIC SURFACE 



-S7 0 1 



-S7 0 2 



-S7 0 3 



± 

RETURN ^ 



Patent provided by Sughrue Mion, PLLC - http://www.sughrue.com 
29 



EP 0 639 753 A2 



FIG. 13 



^ START ^ 



FIND PEAK AND VALLEY FOR Rm( A) AT 200 - 
350nm, CALCULATE SUM (od) OF EFFECTIVE 
OPTICAL PATH LENGTHS OF d3 & d4 
od = 0. 2 5 x (m- 1 ) xw 1 xw2/ (w2 -w 1 ) 



-S210 



CALCULATE MAXIMUM d3max OF d3 

AND MAXIMUM d4max OF d4 
d3max= (1+r) xod/<n3> 
d4max= (1+r) Xod/<n4> 



I 



-S21 1 



CALCULATE FILM THICKNESS PITCHES Ad3, Ad4 
Ad3=0. IX (4x<n3>x (350-250)/ (350x200) ) K 
Ad4=0. lx (4x<n4>x (3 50-2 50) / (3 50x200) ) S212 



± 



DETERMINE Emin 



CHANGE 


d3 


FROM 


0 TO 


d3max 


BY 


Ad3 


: >C 


V 


CHANGE 


d4 


FROM 


0 TO 


d4max 


BY 


Ad4 



hS213 
S214 



S2 i^(l-r)xodg 

<n3>xd3+<n4>xd4 
^ (1+r) xod^ 



YE S 




NO 



Emin«-E 
d3min*-d3 
d 4min*-d 4 



S21S 



S220 



30 



EP 0 639 753 A2 



FIG. 14 




S221 



ASSUME dl=0, d2=0, d3 AND d4 ARE NOT KNOWN, 
IDENTIFY COMBINATION (d3,d4) HAVING MINIMUM 
DEVIATION BY GAUSS NEWTON METHOD 



S222 



CHANG d3 & d4 BY ± 2nm 
FOR NINE CO) 


< .y... 
TO DETERMINE DEVIATIONS 
JBINATIONS 




^S223 


QUADRATIC SURFACE APPROXIMATION USING LEAST 
SQUARE METHOD ON NINE DATA 




^S224 


FIND d3 & d4 HAVING SMALLEST QUADRATIC 
SURFACE, UPDATE d3 & D4 



NO 




S225 



31 



EP 0 639 753 A2 



FIG. 15 



d3 
d 2 

d r 

dO 



S i02 



S i 



S i 02 



S i 



3 
2 
1 
B 



> 1 0 



FIG. 16 



d4 



d3 



d 2 

d r 



dO 



S i Nx 



S i 02 



S i 



S i 02 



S i 



4 
3 
2 
1 
B 



> 1 0 



32 



